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Abstract 



The spacetime singularities play a useful role in gravitational theories by distinguishing physical 
solutions from non-physical ones. The problem, we studying in this paper is: are these singularities 
stable? To answer this question, we have analyzed the general problem of stability of the family of 
QQ . the static spherically symmetric solutions of the standard Einstein-Maxwell model coupled to an ex- 

tra free massless scalar field. We have obtained the equations for the axial and polar perturbations. 
The stability against axial perturbations has been proven. 



X/^ • PACS number(s): 04.20.Cv, 04.20.Jb, 04.40.Nr, 04.80.Cc 



1 Introduction 



QJ , Recently there has been considerable interest in so-called "dilaton fields", i.e. neutral scalar fields 

' whose background values determine the strength of the coupling constants in the effective four- 

dimensional theory. However, although the scalar field naturally arises in theory, its existence from 
^ ', the point of view of the general relativity is quite problematic. It has been shown that including 

^ I a scalar field in the theory leads to a violation of the strong equivalence principle and modification 

of large-scale gravitational phenomena Q. The presence of the scalar field affects the equations of 
motion of the other matter fields as well. Thus, for example, solutions which correspond to a pure 
electromagnetic field appear to be drastically modified by the scalar field. Such solutions were stud- 
ied in ||2|-Q, where it was shown that the scalar field generally destroys the horizons leading to the 
singularities in a scalar curvature on a finite radii. Special attention has been paid to the charged 
dilaton black hole solution ^ . Thus analysis of the perturbations around the extreme charged dilaton 
black hole solution performed in demonstrates the analogy of the behavior of the black holes and 
elementary particles in the sense that there exists an energy gap in the excitation spectrum of the 
black hole. 

From the other side, an interesting way to treat the problem of appearance of the spacetime 
singularities is to develop a theory of gravity including an extra spatial dimensions Q. It turns 
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out that certain singularities can be resolved by simply passing into a higher-dimentional theory of 
gravity for which spacetime is only effectively four-dimensional below some compactification scale ||^. 
Moreover, while studying the decay of magnetic fields in Kaluza-Klein theory was argued in that for 
a physical four-dimensional magnetic field there are two ways it can decay: either by producing single 
naked singularities into which space "collapses," or by producing pairs of monopole-anti-monopole 
pairs which accelerate off to infinity. Since many currently popular unified fields theories include 
an extra spatial dimensions, it is important to ask: could these singularities be stable in our four- 
dimentional world? Although, it was shown that static spherically symmetric solutions for the related 
case of Einstein-Klein-Gordon equations with a quadratic self-interaction term are unstable [Q, it 
would be interesting to study this problem for the general case of the Einstein-Maxwell-scalar system. 

In this paper we consider the problem of stability of the general class of static spherically-symmetric 
solutions of the standard Einstein-Maxwell model with an extra free scalar field (j) with four-dimentional 
action taken to be: 

S = -r^ I dx^^{R - 2g™"V^</.V„0 + F^) , (1) 

where Fab = ^[a^b] is the usual Maxwell field. The geometrical units c = 7 = 1 are used through the 

paper as is the following metric convention (H ). The fields equations corresponding to action 

(1) are easily calculated to be: 

Rmn = '^^m(l)^n4> — '^FmkFn'' + -gmnF'^, (2a) 

g^'VaVbCp = 0, VaF"^ = {2b) 

The general static spherically symmetric solution to system of equations (2) is well known ||2|-||3| 
and it might be given by the following relations 

ds'^ = u{r)dt^ — v{r)dr'^ — w{r)dil., (3a) 



v{r) = —-- = q^{r), w{r) = (r^ - M^)g^(r), (36) 
u[r) 



,M=p,(!:^)Vp,(!i±ii)', (3d) 

( 02 \l/2 

2p± = 1 ± ( 1 + I , (Ao = /iVl-4/c2, (3e) 



where fi,k,Q are the arbitrary constants, the prime denotes the derivative d/dr and the usual notation 
is accepted in (3a) for dQ = dO"^ + sin^ Odip'^. The parameter /i is related to physical mass /xq > and 
charge Q by 
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which saturates the bound \Q\ < ^q. In the extreme limit \Q\ = fiQ, and the solution, independently 
on the scalar field, accepts the familiar form of the extreme Reisner-Nordstrom black hole solution: 
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^{r)=0, R = r±fio. 

In some special cases this solution coincides with well known results which will support our future 
conclusions. This class of the solutions (3) describes the exterior region of the black holes and the 
naked singularities^. It is important to ask: could a distant observer study the objects, located under 
this spurious singularity at r = fi? The answer appears to be no. Indeed, let us imagine that the 
observer will try to test this region using the perturbations of the fields involved. Can the infinite 
energy density (and corresponding singularity in the equation for the perturbations) be an opaque 
boundary for the perturbations, or there is a possibility that the perturbations might penetrate under 
the surface r = ji? To answer this question, one might easily show that r = /i is an effectively infinite 
point in the case of Schwatzchild {k = ±1/2, Q = 0) and Reisner-Nordstrom (k = ±1/2) solutions. In 
the vicinity of this surface there are solutions which propagate in both the "in" and "out" directions 
(~ e^'^''*) Then, for a distant observer, the time of the fall is logarithmically infinite because 
900/911 ~ (r — /x)^ {i.e. we have a horizon). This is true even though the distance to the horizon can 
be traversed in finite proper time. In the general case of the solution (3) with non-zero scalar field 
{k / ±1/2 and \Q\ ^ fj-o), there are no "in" or "out" going waves and one can see that because of 
the relation 500/511 ~ — m)'^ with d < 2 the time of the fall is finite (it has no horizon). From 
the analysis presented in this paper we will see that the perturbations will "stop" at the point r = jjL 
(the singular point of the equation for perturbations) and thus the observer will not be able to see the 
singularity. It is worth to note that, although the energy densities for both scalar and electromagnetic 
fields in solution (3) are infinite when r — > /i, one may show by strightforward calculation that the 
energy (and the mass) of the solution remain finite. This suggests that we may consider a small 
perturbations around the solution (3) and linearize the field equations (2). Moreover, it is reasonable 
to expect that the corresponding energy of perturbations will be small (and therefore also finite) as 
compared to the energy of the solution (3). 

In this paper we will study the general problem of stability of the solution (3) which describes the 
"exterior" region of the black holes and the naked singularities. It is reasonable to note that one would 
not expect general solutions with naked singularities to be stable since the total mass can be negative. 
However, the analysis presented here will show that the solution (3) is stable at least against axial 
perturbations which, in the light of the results of ||5|, [^], makes this reaserch specifically interesting 
for the general case of Einstein-Maxwell-scalar system, superstrings and Kaluza-Klein theories. The 
outline of this paper is as follows: In the next section we will introduce the definitions accepted 
throughout the paper and will obtain the system of equations for axial and polar perturbations. In 
section 3 we will study the problem of stability of the solution (3) against axial perturbations. In the 
following section 4 we will examine the problem of splitting of the obtained 2x2 matrix equation 
into two independent equations. In the final section 5, we will summarize our results and suggest the 
perspectives of the research on the problem of stability of the static spherically symmetric solutions 
of the Einstein-Maxwell-scalar system. We will also discuss some possible experimental consequences 
of presence of the electromagnetic and scalar fields in the motion of the celestial bodies. 

2 The general system of the equations for the fields perturbations 



^In the special case k — ±1/2 result (3) reduces to the Reisner-Nordstrom solution whose properties has been studied 
extensively |^. 
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It is well-known that in the presence of the non-zero background matter field, the perturbations of 
matter and gravitational fields should be studied simultaneously. Otherwise, the equations of motion 
can appear to be inconsistent and, in any case, they can not be applied to the stability problem of the 
solution under consideration. 

In this section we will obtain the general system of the equations for both axial and polar pertur- 
bations for the system of equations (2). In order to simplify the future calculations, let us introduce 
notation for the perturbations of scalar, electromagnetic and gravitational fields. Due to the symme- 
tries of the background field, this can be done at a rather straightforward way ||lO[| . Indeed, as far 
as the background field does not depend on time, we can write the perturbation of any component 
f{t,x) for any field involved as: 

5f{t,x) = ex.p{iujt)Sf{x). 



Because of the spherical symmetry of the background field we, following |11|, will define the spin- 
weighted spherical harmonics by the equations: 

1 

( - ^2 ± S cot 6* =F ^-^d3)sYira{0, if) = 



\i±s + i){i^s)s±iYi^{e,^), 

where 

Then we may expand all the spin-weighted perturbations through these spherical harmonics as follows: 

(i) . The "scalar" perturbations (i.e. the perturbations of the components without the angular 
indices) can be given by: 

5(t){x) = ^ zir,^{r)Yira{0 , (5a) 

Im 

5AQ{x)=Y,hm{r)Yi^{e,^), 5Ai{x)=Y,nim{r)Yi^{e,^), (56) 

Im Im 

5goQ{x) =^aim{r)YimiO,ip), Sgoi{x) = ^ bimir)Yim{0 , if) , 

Im Im 

Sgii{x) = ^cim{r)YimiO,ip). (5c) 

Im 

(ii) . The "vector" perturbations (i.e. the perturbations of the components with one angular index 
only) can be expanded with respect to the spin- weighted spherical harmonics with the spin ±1 
as: 

1 

6go2ix) ± ^—^6go3ix) = 
sm 

= + '^)[dim{r) =F ieim{r)]±iYim{0, ip), (6a) 



Im 



6gi2{x) ± -^—5gi3{x 



sm ( 



E -^(^ + WUr) =F igim{r)]±iYirrXe,^), (66) 

Im 
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sm 

+ l)[oira{r) T isira{r)]±iYi^{e, (6c) 

Im 



(iii). And, finally, the "tensor" perturbations can be expanded |11| with respect to the spin- weighted 
spherical harmonics with the spin ±2, as: 

5922{x) + — ^55f33(^) = 

sm a 

= [^'™(^) - + ^>im{r)]±2Yim{e, if), (7a) 

Im 

1 2i 

Sg22ix) - —^5g33{x) ± ^-^Sg23{x) = 
sm sm d 



= E V - + + 2)[Q„(r) ^ 2iji^{r)]±2Yim{0,v). (76) 
In order to reduce the effective number of the variables we will perform the gauge transformation: 

Xa^Xa + Ca, (8a) 

where the components of the four-vector (,aii, x) are given by the relations: 



^o{t, x) = Y aim{r)Yim{0, ip), 



Im 

^l{t,x)=Y/3lrn{r)Yira{e,^), 
Im 

^2{t,x)±-^^3{t,x) = 

sm 

= Y -K^ + l)[7«m(r) T iSim{r)]±iYim{e, if). {8b) 

Im 

We will impose the same conditions on the coefficients as in [pX|]: 

_ 1. . _ . 

Tim 2^'™' ^Im Jim) 
<^lm = dim - iiO^lm, Plm = flm - I'lm + llm • (9) 

■w 

Furthermore, we will introduce an additional set of convenient notations (with tildas) given by the 
following relations: 

u' 

aim = aim + 2iujaim — Pim — , 

V 

u' 



km = km + i^Plm + a'lm " "im — 



V 



u 
I 



Clm — Qj72 -|- '2(5im Plm > 

V 



u' 

blm = hm + '^i^^aim — Plm. — , 

V 



- ~ , X' X ^ 

9lm — gim + — Olm , 

W 



him — hiyji -\~ Plm ) 
V 



Zlm — hm — ^ {^flm — ^im + ^'™^) ' ^^^^ 

The notations introduced above significantly simplify the future analysis of the perturbations of 
the equations of motion. Thus, by expanding the equations of motion (2) over the field variations and 
then separating the terms with different angular dependence (i.e. terms, proportional to QYim{0,(p), 
±iyim{d,^), ■ ■ ■), one can easily find the correspondent equations for the perturbations. In particular, 
from expressions for the components R22, R23 and given by the equations (2a) we will obtain the 
following relations: 

a'lm = 9lm[ + ii^-eim- (Ha) 

Another equation might be obtained from the expressions for the components R12 and i?i3 (2a), 
namely; 

= e/m— + ~gim^ + 1) - 2 - UJ^-] + 2simA'o. (116) 

W VjJW \ u J 

And finally from the second equation in (26) one may find the last equation: 

^'Im = ^'im ( ~ ) + (11c) 



2v 2u 



-Sim— \ + -U) — + 6;^ - eim l^^gim 

W \ U J \ W J 



U 



Thus, we have obtained three independent components of the perturbations. These components 
are interacting only with each other [10| and hence they have no influence on the other components. 



This is the trivial consequence of the fact that these components are axial, i.e. when the spatial 
coordinates are inverted, their transformation rules appears to be — (— )' rather then (— )'. 
Analogously, the general system for the polar perturbations takes the form: 

u' w' \ ~ u f 1(1 + 1) 



«Zm = «im + ^ + ^l^' ^ ^ ~ 

2u 2w J UjJW \ 2 u 



„ u ~ u I u w ^ - f , 

-ClmTT + "-«m— 7^ 7^ + Hm2AQ - ZimUCj) , 

2v w \2u 2w , 



6 



Clr. 



V 
U 



h'. 



Im 



U 



W 



2u 2w 



A' 



1(1 + 1) 

2iu 



w 



2v 



A' 



~l ^0 --U 

^Ira ~ ^Ira 7; — h him ~ + "f^h 

2u 2v w iLOW 



1(1 + l)-u;^ 



u 



V U 

2v 2u 



+ k 



lUJV 



im' 



U 



w' ^ u' 
w 2u 



aim 

U(f)' 

l{l + l) 



1(1 + 1) 



w 2u 



+ 



3u w 
4uw 



+ 



- 1 



+ 



W(p 

km i 
uw(p' 2uj 



V ^ u / w 
w 2u \2w 



u 
2u 



+ 



1/2 
^0 



— u 



+ 



1(1 + 1) , 



I 2 / \ r ^' 
u - u w 1 + kim -7 + ni,„ . 

W U(p lljJlLKp' 



u u 
l(l + l)A'o 



(12) 



This is the system of independent equations for polar perturbations. All the other equations appear 
to be a consequence of them. 



3 The stability against the axial perturbations 

In this section we will concentrate on the stability of the solution (3) against the axial perturbations. 
In order to approach this problem, we must rewrite the system (11) as an eigenproblem with respect 

to cj^. There exists only one way to combine the first two equations of the system (11) into a single 
equation of the Shroedinger-type where playing the role of energy. To show this, let us define the 
following combination: 

Clm = Meim + ^29lm- 

It is strightforward to check that the equation for C'/^ acquires the form of the Shroedinger equation 
with that substitution only when Ai = 0. 

To present the system of equations for gi^ and s;^ in the hermitian form, it convinient to introduce 
new functions t/ji and ■02 as follows: 

9lm = il^iw^^'^(-) ' Situ, (13a) 



Sim = 4.^2(1) ^Jl(l + 1) - 2. (136) 

(::)"-^'<^>(::)''(nr r,)(::)-' 



Then the equations for the column ^|) = [Jplj takes the form: 
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where the functions a{r) , b{r) , d{r) and weight p{r) are defined form the relations (3) and given as 
fohows: 



„(,) = 5f(l±f)+^?i^. ,14.) 

(r^ — /i^j^ — q[rj 



qyr) (r^ — fi^)-^'^ 



(r^ — /X 



p{r) = q{r)'^. (14e) 

In order to prove the stabihty of the solution (3) with respect to axial perturbations, following 
Wald |12| it is necessary to show that: (i) the spectrum of the differential operator (14) is positive, 
(ii) the differential operator (14) is not only hermitian, but also self-adjoint. 

Concerning the eigenvalues of the operator (14), by straightforward computation one can prove 
that for any p, < r < co both eigenvalues of the potential matrix are positive and then that the 
eigenvalues of the total operator (14) are also positive. 

Now let us analyze the self-adjointness of (14). The boundary conditions at the spatial infinity 
(r +oo) are fixed by the means of the standard procedure: 

dr q{r)'^ip^ip < oo (15) 

and need no further consideration. However, the condition (15) with r ^ fj, permits both possible 
asymptotics for the function ^(r): 

^i(r) = const ■ (r - ^)i/2±(s-i)^ ^^q^-^ 

Mr) = const ■ (r - ^)i/2±(^-i)/2, (166) 

where s = \2k\. By using the condition (15) one might immediately conclude that because of the 
relation: 



q{r fi) = const ■ {r — fi) 

the positive sign in (16a) is forbidden. It means that in order to make the differential operator in (14) 
self-adjoint, we must impose some reasonable boundary condition that will suppress one of possible 
asymptotics of ip2 (166) for r ^ ii. The appropriate restriction appears to be quite natural: to impose 
the condition of finiteness of the energy of electromagnetic perturbations. For the positive sign in the 
condition (166), the energy density is proportional to (r — /u)^^, and the corresponding total energy 
becomes infinite. 

This result completes the proof of the stability of the solution (3) against the axial perturbations 
with the finite value of the initial energy. 
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4 The separation of the equations for V^i(r) and ip2{r) 



Now the question arises whether the system (14) can be spUt into two independent equations for 
some hnear combinations of ipi, ip2 and Let us obtain the general condition on the coefficients 
a, 6, d and p in (14) which wih permit one to say whether 2x2 system can be spht or not. 

The weight p can be eUminated from the equation (14) using the substitution r — > f(r). Thus it 
appears to be sufficient to study the case p = 1 only. With this restriction one will get the following 
equation: 

Let us suppose that there exists set of coefficients r]i, ri2, rjs and 774 that the linear combination 



C = mi^i + V21P2 + mi^'i + ^4V'2 (is) 

satisfies the following equation 

C" + ^'C + J^C = o (19) 

Note that the coefficients in the substitution (18) should not depend on w, otherwise the problem of 
the construction of the coefficients 771, 772, ?73 and r/4 becomes trivial. Moreover, the result obtained in 
this case appears to be practically useless. Indeed, due to the "shadowing" produced by the functions 
^1; '?2j and 774, the behavior of function C, after the substition (18) will not be directly connected 
with the behavior of the initial function ip. 

By comparing the equations (17)-(19), and separating the terms proportional to ■01, ■02, V'i, '02 
and o;^, one can easily find that 



773 = const, r/4 = const. (20) 

It should be noted that the presence of the arbitrary constants 773 and 774 corresponds to the orthogonal 

rotation with the constant coefficients in the {1^1,1^2) plane performed before the definition of the 
function C given by (18). Consequently, keeping in mind the possibility of the preliminary constant 
orthogonal rotation, we can choose 773 = 1, 774 = 0. Then we can find explicit expressions for 771 and 

V2- 



1 b(r) 1 /■ / s , f (b(r) f a(r)dr)dr , 

m = — ^ , ; ; , + - / a(r)dr - ■' ^ ^ , \ , , (21a) 

' 2jb{r)dr 2j^' 2 J b{r)dr ' ^ ' 

m = \j b{r)dr. (216) 



And, finally, the last equation of the system might be presented as: 

air) + 2d{r) + const^ = \-^j^^ + 



^( f h( \^ ^2 , const2 i J [a{r) J b{r)dr]drf 

+^ J ^^'^'^'^ + UKr)drf + 2Ub{r)drf " ^''^ 

This equation is the consistency condition. It means that if it is fulfilled, then the coefficients 771, 772, 
7/3 and 774 given by (20)-(21) satisfy the equations (17)-(19) siniTdtaneously. By restoring the weight p 
(14e), we might obtain from (22) the general form of the consistency condition as follows: 
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a(r) + 2d{r) + consti ■ q'^{r) = (23a) 

a" 1 fa'\'^ 2 2, const2 { J a{r)a{r)drf 
' + p a -\ 5 h 



where the function c7(r) is defined by 

1 f b{r)dr _ 1 /" b{r)dr 
Let us clarify the nature of the equation (23a). It is well-known that any equation of type (14) has 



a "dual" equation |13]. The simplest way to obtain the dual equation is to rewrite the system (14) as 



the 4x4 system of the first-order equations: 

(24) 



i^V _ / M - {q'/q)E iuoq^E \/'ip 



xJ V We -M-{q'/q)Ej \x 

where M is 2 x 2 matrix, E is identity 2x2 matrix. Note that the matrix M always exists because 
it can be directly constructed as 



M = A'A-^, (25a) 



with the matrix A given by 



(256) 



where '0^^-' and ^^^-^ are two linearly independent solutions of the equation (14). The straightforward 
calculation permits us to verify that the function ij: from the equation (24) satisfies the system (14), 
and 

\^ = {M-{q'/q)E)\[M-{q'/q)E)' . (26) 
Similarly, the function x is governed by the "dual" equation 

, . ,2„4 (Xi\ , ( d + d b\ f xi 

where 



d + d b\ X ^\ 2 

6 d 



(m + {q'/q)Ey - (m + {q'/q)E)' . (276) 



One can verify that equation (23a) is equivalent to the condition 6 = 0. It means that equations 
(17) can be separated only if the "dual" system has a diagonal form. 

The straightforward verification of the consistency of the equation (23a) for the values of a, b,d and 
p given by (146) — (14e) leads to rather complicated calculations. From the other side, by analyzing the 
asymptotic behavior of condition (23a) in the limit r — > ^u, we might conclude that (even taking into 
account the preliminary constant ortogonal rotation in (-01,^/^2) plane) equation (23a) can't be fulfilled 
in the limit r — > /i. It leads us to the conclusion that the system (14) is "essentially" two-dimensional 
and can not be split into two independent equations. 
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5 Discussion. 



We have analyzed the problem of stability of the exact solution of the standard Einstein-Maxwell 
gravity coupled to an extra free massless scalar field. It was shown that, although the solution (3) 
contains naked singularities, it is stable at least against axial perturbations. The problem of the 
stability of this solution against polar perturbations is much harder to analyze. One unexpected 
complication of these studies is that the differential operator for corresponding 3x3 eigenproblem 
appears to be non-Hermitian. However, this research is currently in progress and the obtained results 
will be reported in a subsequent publication. 

Anticipating the possible questions we would like to note that the correspondence of our analysis to 
the existing results on the perturbations of the Reisner-Nordstrom solution is not quite strightforward. 
The reason for this comes from the conclusion that no smooth limit of our analysis exists for 
and, although we can reduce our relations to the case k = ±1/2, the final results will be degenerate at 
the point r = fi. To show this, one may examine the case 500/511 ~ — A*)'^ with d = 2 and see that 
the frequency term in the equation (14a) is infiuencing the asymptotic behaviour of the function iplr) 
at the vicinity of the surface r = ^ Moreover, the matrix structure of the eigenproblem becomes: 
E ■ a(r) + (T3 • b{r) ■ consti + o"i • 6(r) • const2, where &i and are corresponding Dirac matricies. 
As a result, the equations for perturbations in the case of the Reisner-Nordstrom solution appear to 
be split into the following three groups (reconstructing the already known result [^): (i) the scalar 
perturbations, which don't interact with other perturbations; (ii) the two independent modes of the 
axial perturbations, where both the gravity and the Maxwell field are mixed together, and (iii) the 
two independent modes of the polar perturbations, which also contain a mixture of gravitational and 
electromagnetic fields. 

Concluding this part, we would like to note that the stability of the general solution which contains 
the naked singularities is well fit to the scenario proposed in the multi-dimentional extensions of the 
general relativity. Thus, under the certain circumstances, the Kaluza-Klein vacuum may decay by 
endlessly producing naked singularities. This process from the five-dimensional point of view corre- 
sponds to Witten's "bubbles of nothing" which must eventually collide [0], and so in four-dimensions 
the singularities will coalesce. However, it should be emphasized that we have explored just the four- 
dimensional solutions and the further analysis of this problem should include the non-trivial coupling 
of the scalar field in the higher dimensions. 

The proven stability of the exterior static solution (3), makes it interesting to study whether the 
scalar and electromagnetic fields might be detected through the space gravitational experiments. Thus, 
following the standard procedure of the PPN formalism |jl^, one will find that only the parameter (5 
deviates from its general relativistic value, namely 

where is the Newtonian mass of the source and we have restored the dimensional constants 7 and 
c. This result coincide with one for the Reissner-Nordstrom solution and as long as the metric (3) in 
post-Newtonian limit doesn't contains parameter /c, the scalar field (defined by the action (1)), can 
not be detected from the data processing of the modern relativistic celestial mechanical experiments. 

The second term in the expression (28) represents the ratio of the electrostatic energy contribution 
in the gravitational field produced by the same charged massive body. The presence of this term 
might lead to an observable discrepancy in the motion of the celestial bodies. For example, it will 
contribute to the Nordtvedt effect, which was extensively studied in the Moon's motion [p!4[. In the 



^For the general case with non-zero scalar field {k ^ ±1/2 and \0\ 7^ /io), as we saw, d < 2 and this influence is 
absent. 
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recent analysis of data obtained in Lunar Laser Ranging which was carried out to detect the Nordtvedt 
effect, a very tight ]l5[ | hmitation on the parameterized post-Newtonian parameter (3 was obtained: 



P = 0.9999 ± 0.0006 (29) 

This result suggests that within the accuracy one part in ten thousand, the contribution of the ratios 
of electrostatic to self-graviational energies presented by (29) for both Moon and Earth is negligable 
small. 

Since gravity attracts positive and negative charges equally, then the matter accreted on a massive 
astrophysical object will be nearly neutral. For the case of the celestial bodies with the gravitational 
energy dominating over the electromagnetic one, the parameter /3 might be presented as: 

where the mass of the star /xq was expressed in terms of the solar masses Mq: /iq = uMq. The con- 
straints imposed on new weak forces from the behavior of the astrophysical objects gives for the maxi- 
mum possible electric charge Qmax carried by celestial bodies the following estimation: Qmax 



|16|. This gives the following estimation for the electrostatic energy contribution in the parametrized 



post-Newtonian parameter /3: 



A/3 = /3- 1 < ^ X 10"^ (31) 



Unfortunately, even with n = 1 this result gives practically unmeasureable value for the contribution 
of the electrostatic energy of the charged astrophysical body to the generated gravitational field. 
According to this result, the detection of the electrostatic field contribution in the relativistic celestial 
mechanics experiments performed in the weak gravitational field is presently impossible. 

Thus, we have shown that the infiuence of both the electromagnetic and the scalar fields (given 
by the action (1)) on the motion of the asprophysical bodies is practically unmeasureable in modern 
gravitational experiments. However, a wide class of multi-dimensional theories of gravity with an 
arbitrary number of massless scalar fields coupled to the usual tensor gravitational field has been 
recently considered in [^. This investigation was performed in order to analyze the cosmological 
consequences of an inclusion of the multi-scalar field terms in the theory. As a result, the authors of 
this paper illustrated that although these theories might have coinciding post-Newtonian limits with 
general relativity, they predict non-Einsteinian behavior of the stellar objects in a strong gravitational 
field. In particular, it was noted that this discrepancy will lead to observable effects, for example, for 
the binary pulsars. This result makes it specifically interesting to study the multi-scalar field extensions 
of the general Einstein-Maxwell-scalar model together with the condition to meet the experimental 
constraints based on the tests of general relativity performed to date. 
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